function f = ind_shares(expdelta_cf,expmu_cf)

ns=size(expmu_cf,2);
expon=repmat(expdelta_cf,1,ns).*expmu_cf; 
denomi=repmat((sum(expon,1)+1),length(expdelta_cf),1);      

f=expon./denomi;
